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Abstract 

For the free models of statistical mechanics on torus, exact asymptotic ex- 
pansions of the free energy, the internal energy and the specific heat in the 
vicinity of the critical point are found. It is shown that there is direct relation 
between the terms of the expansion and the Kronecker's double series. The 
latter can be expressed in terms of the elliptic 0-functions in all orders of the 
asymptotic expansion. 
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I. INTRODUCTION 



It is well known that the singularities in thermodynamic functions associated with a 
critical point occur only in the thermodynamic limit when dimension L of the system under 
consideration tends to infinity. In such a limit, the critical fluctuations are correlated over 
a distance of the order of correlation length ^buik that may be defined as the length scale 
governing the exponential decay of correlation functions. Besides these two fundamental 
lengths, L and ^buik, there is also the microscopic length of interactions a. Thermodynamic 
quantities thus may in principle depend on the dimensionless ratios Cbuik/-^ and a/L. The 
Finite-Size Scaling (FSS) hypothesis [0 assumes that in the scaling interval, for temperatures 
so close to the critical point that a <^ ^buik ~ L, the microscopic length drops out and the 
behaviour of any thermodynamic quantity can be described in terms of the universal scaling 
function of the scaling variable t = C,huik/L. However, non-universal corrections to FSS do 
exist. These sometimes can be viewed as asymptotic series in powers of a/L. 

Two-dimensional models of statistical mechanics have long served as a proving ground 
in attempts to understand critical behavior and to test the general ideas of FSS. Very few 
of them have been solved exactly . Notably, the so-called free models that can be treated 
in terms of non-interacting quasi-particles on the lattice. The Ising model (which is 
equivalent to the model of free lattice fermions) and the Gaussian model (which can be 
viewed as the model of free lattice bosons) are the most prominent examples. 

For all free models, exact asymptotic expansion of the free energy on an infinite cylinder 
of circumference L can easily be obtained by direct application of the Euler-Maclaurin 
summation formula 0,^. However, derivation of such an expansion on a torus of area S 
and aspect ratio p is much more difficult problem. For the Ising model on torus, such an 
expansion has first been studied by Ferdinand and Fisher 0. Starting with the explicit 
expression for the partition function 0], they have calculated two leading terms, /bulk and 
/o(p), of the expansion 

oo 

Ft=tc{p, S) = /bulk S + /o(p) + E fpip) (1) 

p=i 

Their calculations have recently been pushed forward to get two next sub-leading terms 

Eg- 

It is the purpose of this paper to derive all terms of the exact asymptotic expansion of 
the logarithm of the partition function on torus for a class of free exactly solvable models of 
statistical mechanics. Our approach is based on an intimate relation between the terms of 
the asymptotic expansion and the so-called Kronecker's double series. Besides the aesthetic 
appeal of the exact expansion, there is also physical motivation to study non-universal 
corrections to FSS. The problem is that in numerical simulations of lattice models one usually 
studies relatively small lattices. Therefore, to compare the results of high precision numerical 
simulations to the theoretical predictions one cannot neglect sub-leading corrections to FSS 
||1CI||. Non- universal terms in the asymptotic expansion also provide important information 



about the structure of irrelevant operators in conformal field theory |11 
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II. FREE MODELS OF STATISTICAL MECHANICS 



In this section, we formulate three basic models of statistical mechanics: Ising model, 
dimer model and Gaussian model. These models are often refereed to as "free" since they 
were shown to be equivalent to free fermions or free bosons on the lattice. Partition functions 
of all these models can be written in terms of the only object — the partition function with 
twisted boundary conditions Za^p{ii). Exact asymptotic expansion of the later will be our 
main objective in the subsequent section. 



A. Ising model 



The Ising model is usually formulated as follows. Consider a planar square lattice of size 
M X N with periodic boundary conditions, i.e. torus. To each site (m, n) of the torus a spin 
variable is ascribed, Smn, with two possible values: +1 or —1. Two nearest neighbor spins, 
say Smn and Smn+i contribute a term — J Smn Smn+i to the Hamiltonian, where J is some 
fixed energy. Therefore, the Hamiltonian is simply the sum of all such terms, one for each 
edge of the lattice 



N-l M-1 
n=0 m=0 



(2) 



The partition function of the Ising model is given by the sum over all spin configurations on 
the lattice 

It is convenient to set up another parameterizations of the interaction constant J in terms 
of the mass variable /i = In ^/sh2J. Critical point corresponds to the massless case /x = 0. 

An explicit expression for the partition function of the Ising model on M x torus, 
which was given originally by Kaufmann |12|, can be written as 

^isi„g(/i) = I (v^e^)*'"^ {^|,|(/^) + ^o,i(/^) + Zi^M + ZoM} (3) 
where we have introduced the partition function with twisted boundary conditions 

N-l M-1 

ziM = n n 4 [ sin^ + sin^ + 2 sh^ ; 

n=0 m=0 

Here a = corresponds to the periodic boundary conditions for the underlying free fermion 
in the iV-direction while « = | stands for anti-periodic boundary conditions. Similarly j3 
controls boundary conditions in M-direction. With the help of the identity []13| 



M-1 



4|sh(Moj + i7r/3)|^ = 4 sh^Moj + sin^ vr/? 



n 4[sh2u; + sin2 



m=0 



M J 



the partition function with twisted boundary conditions Z^^p can be transformed into simpler 
form 
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N-l 



n=0 

where lattice dispersion relation has appeared 



Z^M = n 2 I sh [Mu, (^^) + t7T(3] I (4) 



i^tj.{k) = arcsinhy sin k + 2sh. fi (5) 

This is nothing but the functional relation between energy u;^ and momentum /c of a free 
quasi-particle on the planar square lattice. 



B. Dimer model 

A "dimer" is a two-atom molecule. When drawn on a lattice it covers two adjacent sites 
of the lattice and the bond that joins them. The "dimer problem" is to determine the number 
of ways of covering of a given lattice with dimers, so that all sites are occupied and no two 
dimers overlap. If we consider planar square lattice of size 2M x 2N wrapped on a torus 
then the number of dimers must be 2MN and the number of distinct dense coverings of the 
lattice (the partition function) has been calculated by Kasteleyn and Fisher |jl^,|l5|. This 
can be expressed in terms of the same partition function with twisted boundary conditions 
Eq.(|D as 



^ Dimer 



I [Zl ,(0) + Zlm + Zl „(0) - Zl,{Q)} (6) 

Z ^ 2'2 '2 2' ' J 

First two leading terms of the asymptotic expansion of this partition function has been 



obtained by Ferdinand |T6|. Let us also mention that dimers are always in the critical point 



and have no phase transition. 



C. Gaussian model 

Let us now turn to a boson analog of Ising model, which is often referred to as Gaussian 
model. Again, we consider planar square lattice of size N x M wrapped on a torus. To 
each site (m, n) of the lattice we assign a continuous variable Xmn- The Hamiltonian of the 
model is 

N-l M-1 
n=0 m=0 

The partition function of the model can be written as 



ZiJ) = [ e-^(^-) daix) 

J-dMN 



If the measure da{x) in the phase space R*''^^ is chosen to be Gaussian 

Af-l M-l 



n=0 m=0 
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the integration can be done explicitly and the partition function of the free boson model can 
be written in terms of the partition function with twisted boundary conditions Eq.(^ and 
parameterization = 4 ch^ /i as 



^ Gauss 



(/i) = (v^ch/i)*''^ \ZoM] ' (8) 



This model exhibit phase transition at the point /Zc = where the partition function is 
divergent. This is due to the presence of so-called zero mode, i.e. due to the symmetry 
transformation Xmn Xmn + const, which leave the Hamiltonian (|^ invariant. Correlation 
functions of disorder operator in this model have been studied by Sato, Miwa and Jimbo 

The reason why this model is often considered as boson analog of the Ising model is that, 
one can choose another measure in the phase space, which makes this model equivalent to 
the Ising model considered above 



N-l M-1 



dai3ing(x) = 2-*^^ n n ['^(^-- - 1) + ^i^rnn + 1)] dx. 



n=0 m=0 



where 5's are Dirac (5-functions. With such a definition the variables Xmn can actually take 
only two values: +1 or —1, so that a;^„ = 1. In this case integration can be replaced by 
summation over discrete values of Xmn = ±1 and the Hamiltonian ([^) coincides with the 
Hamiltonian of the Ising model (0) up to a constant. 



III. ASYMPTOTIC EXPANSION IN TERMS OF KRONECKER'S DOUBLE 

SERIES 

In the previous section it was shown that partition functions of three basic free models 
of statistical mechanics can be expressed in terms of the only object: the partition function 
with twisted boundary condition Za^/^i^). In this section we shall obtain exact asymptotic 
l/5'-expansion of the partition function near the critical point. For reader's convenience, 
all the technical details of our calculations and the definitions of the special functions are 
summarized in the appendices attached to the paper. 

First of all, let us mention the symmetry properties of the partition function 
From its definition (|^) one can easily verify that it is even and periodic with respect to its 
arguments a and /3 

ZaAf^) = Za-(}{fJ') = Z^aAl^) 
Za,l3{fJ') = Zi^a,f3{fJ') = Za,l+p{fi) 

These imply that twist angles a and (3 can be taken from the interval [0, 1]. Then, one can 
note that for all twists {a,P) ^ (0,0) the partition function Za,i3{fi) is even with respect to 
its mass argument fi. Hence, near the critical point (/i = 0) we have 

Za,M = ZaAO) + C Z'^AO) + ■■■ («, /3) 7^ (0, 0) (9) 
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The only exception is the point where both a and (3 are equal to zero. This case has to be 
treated separately since at this point the partition function turns to zero. As a result, we 
have 

Z,M = /i 4,o(0) + ^ Z^^(O) + . . . (a, (3) = (0, 0) (10) 
In what follows notation Za^p^jj) will imply (a,/?) ^ (0,0). 



A. Asymptotic Expansion of Z^^pifi) 



Considering the logarithm of the partition function with twisted boundary conditions, 
Eq.(HD, we note, that it can be transformed as 



N-l 



N-l 



N 



lnZ,,^(0) = M5]a;o 

n=0 ra=0 

The second sum here vanishes in the formal limit M 



-2[a/..o(^^)+*^/3] 



(11) 



oo when the torus turns into 
infinitely long cylinder of circumference N . Therefore, the first sum gives the logarithm of 
the partition function with twisted angle a on that cylinder. Its asymptotic expansion can 
be found with the help of the Euler-Maclaurin summation formula (Appendix ^ 



P \ T3a 



{2p)\ 2p + 2 



;i2) 



Where /Ja;o(x)dx = 27, 7 = 0.915965... is Catalan's constant and are so-called Bernoulli 
polynomials. We have also used the symmetry property, u!o{k) = uq^tt — k), of the lattice 
dispersion relation (^) and its Taylor expansion 



uo{k) = k\ X + J2 



A 



2p 



(13) 



where A = 1, A2 = —2/3, A4 = 4, etc. In what follows, we shall not use the special values of 
these coefficients assuming the possibility for generalizations. 
The second term in Eq.(|ll]) we may transform as 

N—l 



n=0 



00 1 i'W/21-1 
^ — 1 THj I ^ — n 



- Re 



m=l 



n=0 

N-[N/2]-l 



+ 



E 

n=0 



-2m Mojo 



' 7r(ri + l — a) ^ 



-ITTf) 



(14) 



The argument of the first exponent can be expanded in powers of 1/5 if we replace the 
lattice dispersion relation ljo{x) with its Taylor expansion (p!3D 



exp < —2Tcm[ Xp{n + a) + i(3] — 2Tcmp ^ 



A2p / ttV' 



{n + a) 



2p+l 
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Taking into account the relation between moments and cumulants (Appendix 0), we obtain 
asymptotic expansion of the first exponent itself in powers of 1/S 

^-2m[MoJo ( '^'"^"^ )+i7r/3 ] ^ ^-2nm[\p{n+Q)+il3] 



The differential operators that have appeared here can be expressed via coefficients A 
of the expansion of the lattice dispersion relation as (|13| 

A2 = A2 

2 ^ 

A4 = A4 + 3A2 TTT 

aA 

cP' 
Ag = Ae + I5A4A2 — + 15Af 



2p 



Ar 



EE 



r=l 



pi 



Qk 



ki\...kr\ 



Here summation is over all positive numbers {fci . . . /c^} and different positive numbers 
{pi, . . . such that piki + . . . + p^K = p and k = ki + . . . + kr — I. 

The expansion for the second exponent in Eq.([T^) can be obtained along the same lines 
by substitution: a ^ 1 — a. Plugging the expansion of both of the exponents back into 
Eq.(p!4D we obtain 



N-l 



Eln 



n=0 



^ _ ^-2[M..o{^^)+^^p] 



00 . r[iV/2]-l N-lN/2]-l 
Re^—l ^ ^-27Tm[\p{n+a)+i/3\ _^ ^ g-27rr7i[ Ap(n+l-a)+j/3 J 



m=l 



m 



n=0 



n=0 



p=l 



f[N/2]- 



i2p)\ 



m=l I n=0 

N-[N /2]-l 

+ Y _^ 1 _ Q,)2p+1 g-27rm[Ap(n+l-a)+j/3] 

n=0 

In all these series, summation over n can be extended to infinity. The resulting errors are 
exponentially small and do not affect our asymptotic expansion in any finite power of 1/S. 

The key point of our analysis is the observation that all the series that have appeared in 
such an expansion can be obtained by resummation of either elliptic theta function, Oa^T) 
(Appendix [C]), or Kronecker's double series, Kp'^(r) (Appendix 0). Namely, with the help 
of the identities flCW) and fpll) we obtain 



N-l 

Eln 

n=0 



_ 2[Muo{--^)+inp] 



In 



r]{iXp) 



+ 7rApB^ 



pV A,^ Re K^,f_,,{zXp) - B^^, 



{2p)\ 



2p 



(15) 
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Substituting Eqs.([T2D and ([T5|) into Eq.(|Tl|) we finally obtain exact asymptotic expansion 
of the logarithm of the partition function with twisted boundary conditions in terms of the 
Kronecker's double series 



lnZ,,^(0) 



— / ido{x) dx + In 

TT Jo 



TjiiXp) 

^ ~ /vrVV A2p Re K^lS^p) 



(16) 



Note, that Bernoulli polynomials have finally dropped out from the asymptotic expansion 
on torus. This actually means that Kronecker's double series can be considered as elliptic 
generalizations of Bernoulli polynomials. 



B. Asymptotic Expansion of Z'qq{{)) 



As it has already been mentioned, we have to treat the case = (0,0) separately. 

Taking the derivative of Eq.(^) with respect to mass variable p and then considering limit 
/X — >• we obtain 



N-l 



Z^,o(0) = 2V2M n 2 I sh [M^o(f ) 



n=l 



Asymptotic expansion of this expression can be found along the same lines as above. In 
terms of the Kronecker's double series, the expansion can be written as 



In 4,0(0) 



S 



IT Jo 



uJo{x) dx + - In 8pS + 2 In | rj{i\p) 



^ A2p ReK 



0,0 

2p+2 



[iXp) 



{2p)\ 2p + 2 



(17) 



C. Asymptotic Expansion of Z^'^(O) 



The analysis of the Z^'^(O) is a little more involved. Taking the second derivative of 
Eq.(H[) with respect to mass variable p and then considering limit /i ^ we obtain 



N-1 



ReM cth [Mcuo(- 



jn+a) 



n=0 
N-1 



N 

N-1 00 



(18) 



= M ^ uj'^{^^) + 2 ReM J2 E ^0 (^^) e-2-K-o(T)+H 

n=0 n=0 m=l 

where uJq{x) is the second derivative of uj^ix) with respect to p at criticality 

2 

^0 W 



smx + sin^ X 
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Using Taylor's theorem, the asymptotic expansion of the u!o{x) can be written in the follow- 
ing form 



X 



1 + E 



2p 



where k = 2, K2 = —2/3, K4 = 172/15, etc. Again, in what follows, we shall not use the 
special values of these coefficients assuming the possibility for generalizations. 
The first sum in Eq. ( |T8| ) we may transform as 



7V-1 

ME^o(^)=ME/(- 

n=0 



N-l 



(n+g) 
N 



n=0 



+ -E 



+ 



1 



(19) 



where we have introduce the function f{x) = uJq{x) — k/x — k,/{tt — x). This function and 
all its derivatives are integrable over the interval (0, tt). Thus, for the first term in Eq. (|19D 
we may use again the Euler-Maclaurin summation formula (Appendix and after a little 
algebra we obtain 



N—l C pTi 00 

= - / /(x)dx-vrp/.E 

n=0 ^ p=l 



2 \ P-1 

S J p{2p)\ 



00 „\ J^~J^ ^ "Ra 



2p 



(20) 



where J^fix) dx = 21n2 — 41n7r. The second sum in Eq. (|I9|) can be written in terms of 
the digamma function iIj{x). 



N-l 

E 

n=0 



+ 

n + a n + I — a 



[iIj{N + a) + ip{N + 1 - a) - ip{a) - ^{l - a)] 



(21) 



The asymptotic expansion of the digamma function ip[x) is given by (see Appendix 0) 

00 ga 1 

i^{N + a)=\nN- E(-l)^^^ (22) 



Using the symmetry properties of the Bernoulli polynomials B", namely, B 
B2P+1 = ~B2p+i, the Eq. (0) can be rewritten as 



2p 



Bgp and 



Af-l 

E 

n=0 



n + a n + \ — a 



2\nN -Y, 



00 "DO -I 
^2p -L 



ip{Q.) — ?/'(l — a) 



Plugging Eqs. (^0|) and (|23|) back in Eq. ([l^) we have finally obtain 

ip{a) + ipi'^ — «) 



(23) 



N-l 
n=0 



X 2kS 


-/ 


' vr 


2k Jo 



f{x) dx + lnA^ 



2p 



p(2p)! 



(24) 
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Let us now consider the second sum in Eq. (pISl). Note that function uJq{x) can be 
represent as 

4'W^^xp{gg^.*} (25) 

where coefficients e2p and K2p are related to each other through relation between moments 
and cumulants (Appendix jB]). Following along the same lines as in the section (3.1), the 
second sum in Eq. ([T8|) can be written as 

N-l oo 



2 ReM ^ ^ ^o(^) e-2"[*^-°(^)+H = 

n=0 m=l 

= ^^'^ Rc V V [ ^-27Tm[p(n+a)+ip] _^ ^ ^-2nm[p(n+l-a)+il3]\ 

^0,^1 ^^ + « n + l-a J 



oo oo ^ r T r i 

-27rm[ Ap(n+l-a)+j/3j 



71=0 m=l 



(n + 1 — a)e~ 



oo n /^2^\?'-J^ oo ^ R« /^2^\J'-^ 



i-^2p — £2p + ^2p4x as 



The differential operators f22p that have appeared here can be expressed via coefficients 

Q C 

^2 = ^2 



= Li;4 + ScUn 



Let us introduce the function Ra^p{p) 

R.AP) = -*M±|(i^ (27) 

oo OO , 1 1 

+ Re V V i g-27rm[p(n+a)+i/3j _^ ^_27rm[p(n+l-a)+i/3j 

^0,^1 ^^ + « n+l-a 
which ffist derivative with respect to p is given by (see Appendix |y, Eq. (|C^)) 

^i?«,^(p) = -2^1n|^„,^(zp)| + i- (^^ln|^„,^(zp)|^ (28) 

For the cases (a, /5) = (0, 1/2), (1/2, 0), (1/2, 1/2) the second term in above equation is equal 
to zero, and for _Rq^^(p) we obtain 

RaAp) = -2\n\e^Aip)\+CE + 2\n2 (29) 

where C^; is the Euler constant. 

With the help of the identity ( p3D the Eq. (^) can be rewritten as 
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N—1 OO r / r _i_ ^\ 1 

2 ReM ^ ^ ^o(^^) e-2"[*^-°(T)+H 

n=0 m=l 



TT 



'?/'(«) + — a) 



K 



d 



+ 2 + 



,9^ 



In 



ri{ip) 



p=2 



p{2p)\ \ S 



1 P(2p)! I S ) 



(30) 



Substituting Eqs. 
of the Z;^(0) 



^) and (^OD into Eq. (|T^) we finally obtain exact asymptotic expansion 
Zlp{Q) _ 2kS 



TT 



1 

2k Jo 



f{x)dx + \nN + R^^pip) 



^2p ( P 



^,pi2p)\ V S 



r]{ip) 

^2p 



(31) 



Re K-f{z\p) 



D. Asymptotic Expansion of Zqq(O) 
Let us now consider the case (a,/3) = (0,0). Eq. (|10D imphes immediately that 

hm — —— = — ,' , , 
M-o ZoAp) ^0,0 (0) 



(32) 



Taking the third derivative of Eq.(^) with respect to mass variable p and then considering 
limit /i ^ we obtain 



zuo) 



N-l 



2M2 - 1 + 3M ^ u;'^[f) cth [Mu;o{f 



n=l 



(33) 



Asymptotic expansion of the Zqq(O) can be found along the same lines as above. In 
terms of the Kronecker's double series, the expansion can be written as 



^o,o(0) vr 



1 

2k Jo 
d 



f{x) dx + In + Ce — In ri{ip) 



+ 3k K,2PTr- + A2P TTij \nri{ip) - 1 



dp 



dp^ 



K2p°(^Ap) 



(34) 



Expansions (p!6|), (|T^, ( pT] ) and (p^ are the main results of the paper. Kronecker's 
double series K^'^ with a and (3 taking values and 1/2 can all be expressed in terms of the 
elliptic ^-functions only (Appendix |^) . 
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IV. ASYMPTOTIC EXPANSION OF THE FREE ENERGY, 
THE INTERNAL ENERGY AND THE SPECIFIC HEAT 



After reaching this point, one can easily write down all the terms of the exact asymptotic 
expansion (P of the free energy, F = — In Z, at the critical point for all three models under 
consideration. For the Ising model and the Gaussian model we have found that the exact 
asymptotic expansion of the internal energy, U = — din Z{J)/d J, and the specific heat, 
C = d'^ln Z{J)/dJ'^, at the critical point can be written in the following form 

oo 
p=0 

oo oo 

c = cbulk S + C2p(p) s-^'+-^ + J2 C2p+i(p) s-p 

p=0 p=0 



1. The asymptotic expansion of the free energy of the Ising model is 

r- 27 
/bulk = - In V2 



TT 

Mp) = - In- 



flip) 

/2(P) 



27] 

n'p' 1(91 + 91 + 91) + 9,9,9, [9l9l - 9^1 - 9l9l] 



180 92 + 93 + 9, 

vrV^ 929;9, {91(91 - 9tY + 9l{9l + 9lf + 9l{9t + 9tf} 
18432 {92 + 9, + 9,Y 

Slvr^^ 9t{9l - 91) + 9\{9\ - 9%) + 9\{9l - 9%) \ 

92 + 9, + 9, + ^% 

vrV 029,9, [9l{9l - 9l) + 9l{9l - gj) + 9l{9l - 9l)] ( ^ . J \ 

92 + 9, + 9, 

3l7rV5 et9l [Wi9l{9, + 9,) - 9l{29l + 9l - 
72576 92 + 9, + 9, 

vrV 9l9t [929M - 91) + 929,(91 - 9\) + 29\9,9,{9l - 9\) - ^19"^^ 
4536 92 + 9, + 9, 



Wbulk - 

^o(p) = -2a/p 



u\{p) 

Mp) 



V2 

^2 ^3 ^4 

^2 + ^3 + 9, 

7rV/2 ^^2^304(^2' + ^3 + 



48 {92 + 9, + 9,)^ 
vrV^/^ 919191 [9l{9t - 9i)' + 9l{9t + 9^)^ + 9l{9\ + 9 



4608 {92 + 9, + 9,Y 



12 



^ ttV^/' ^2^3^4 [23(^f + + ) + 8^1^1(5^9 + 591 - 8^1 + 20j0l + 20l0i)] 



\2 



9216 {02 + 03 + 0i) 

^'p'" 02030, [0m - et) + em + 01) - 0m + ^p] (, ^ ^ 
-m {02+0s+0.r V ^ ^ % 



Cbulk(p) 



In 



TT 



/ - + In + Cij - - 



-4p 



0203^4 



16 E 



f=2 1^ 



^02 + 03 + 04] n 02 + 03 + 04 



Clip) 



ttV (^i - ^3)^2^3 In I + {01 - 01)0403 In I + {01 - 0l)0,02 In | 

6 (00. + 03 + 0aY 



^ tt'p^ 0l0t{202 - 03 - 04) 

02 + 03 + ^4 
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+ 



np 01 



el 



^3 

e3{el - 



TT^p^ elelem + et + el) 

12 {02 + 03 + 04)^ 
2020a{01 -01) . d 



02 + 03 + ^4 

, ,T^0_2Mm±3±3i 

24^2 (^2 + ^3 + ^4)2 



1 + 4p— ln^2 



We have also used the following relations between derivatives of the eUiptic functions 

I- In ^3 = 7^4 + I- In ^2 and |- In ^4 = j^l + |- In ^2 
op A op op A dp 

Note, that with the help of the identities 



^ In ^2 = -\elE, and = "^elel - lelE 



one can express all derivatives of the elliptic functions in terms of the elliptic functions 
^2, ^3, ^4 and the eUiptic integral of the second kind E 



2. Similar expansion of the free energy of the dimer model is 

/bulk = 

TT 

/o(p) = -ln'^ 



.1 

ei + 0i + 0i 



hip) 
Up) 



27^2 

ttV^ Uof + 0f + 0f) + 010101 \eiel - eiei - 010^ 



90 01 + 01 + 01 

ttv ei (101 + 0i0t) + 01 li0i - 0i0ty + 01 (101 - 0i0iy 



16200 01 + 01 + 01 
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+ 



+ 



+ 



16200 



756 

SIttV 
12096 



DlO I nW I qW\ I olcaca \q2q2 n2n2 



+ 



el + el + el 



01) 



+ ei{ei-ei)] + U'2ef 



Of) 



01 + 01 + dl 



Of[ 



t - 01 



+ Ofi 



n+ot)-of{ 



n + ot) 



01 + 01 



756" 



oloi m 



n - 01) 



I + 01) 



1 + In ^2 



om+oi)] 



01 + + el 



1 + 4p^ In ^2 



3. Finally, the asymptotic expansion of the free energy, the internal energy and the specific 
heat of the Gaussian model after subtraction of the zero modes. In fi\/8S, \/2S/ p and 
S{1/ fj, + l//i^), respectively, can be written as 



/bulk = In a/2 

/o(p) = In Vpr?^ 




For internal energy there are no finite-size correction terms [ubuik 
p). 



0, Up = for all 



, , 8 /, , ^ , ^ 

Cbuik(p) = - In W - + In h - - + In 

TT y y P TT 4 y 

Co(p) = 

2 / 9 92 \ 

ci(p) = 3 ( l + 4p— In r/ + 4p2— In r/j 

C2(p) = 



V. SUMMARY 



In this paper, we have derived exact asymptotic expansion of the partition function with 
twisted boundary conditions at the critical point, Eqs. (0) and (|T^. As an application 
of this result, we have obtained exact asymptotic expansion of the free energy, the internal 
energy and the specific heat for a class of free exactly solvable models of statistical mechanics 
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on torus. Moreover, the partition function of the dimer model on the Klein bottle can also 
be expressed via the partition function with twisted boundary conditions, namely .^Kiein = 
(0) |]18|. Exact asymptotic expansion of the latter can immediately be written down 



Zi 1 

4 



with the help of our general formula (0). An interesting question is: whether this is also 
the case for other free models? This, however, is the problem for the future. 



APPENDIX A: EULER-MACLAURIN SUMMATION FORMULA. 

Suppose that F{x) together with its first 2m derivatives is continuos within the interval 



[a,b). Then the general Euler-Maclaurin summation formula states [|T9l 

n=0 k=l 



^ F{a + nh + ah) = - [ F{t) dr + ^ -rrBfc(a) {F^'-'\b) - F'^''-'\a)) - Rp{a) 

(Al) 

where p < 2m, 0<a<l, h = {b — a)/N and reminder term Iip{a) is given by 

hP ri , r^-i 1 

Rp(a) = — / Bp(a-r) <^ ^ F(f)(a + n/i + r/i) dr (A2) 

P- -'^ [n=0 J 

Bp(a) are so called periodic Bernoulli functions which are defined as follows 

B„{a) = \ V (A3) 

These functions have singularities at the integer values of a and on the interval a G [0, 1] 
they coincide with so-called Bernoulli polynomials 

Bi(a) = a-\ 
B2(a) = — a + I 
B3(a) = - + 
B4(a) = - 2a=^ + - ^ 

The generating function for the Bernoulli polynomials is 



^ - = 1 + y -r B„ 



a) 



Fourier transform of the generating function gives the important identity 

p2mza oo -r +00 —2-Kina 

e2-- - 1 27,1 ^t^^ z + n 

It is well known that Euler-Maclaurin summation formula is closely related to (generally 
speaking divergent) asymptotic series. For further discussion of the properties of these 
series the interested reader is referred to the book of Hardy [pO| . 
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In this paper we are mainly interested in sums of the form 



N ^' \ N 



The asymptotic expansion of the sum ( [A5|) in the hmit N oo can be obtained from Eq. 
( |Al| ) by setting (a = 0,6 = tt). If we assume that all the derivatives of f{x) are integrable 
over the interval (0,7r), i.e., the integral in Eq. is finite, we can formally extend the 

sum in Eq. ([Al|) to = cxd and drop the reminder term Rp(a). In this case we can write 



the asymptotic expansion of the sum as follows 



n=0 \ / " •''^ " k=l 

Note, that we also use abbreviated notation for Bernoulli polynomials: Bp (a) = B^. 

APPENDIX B: RELATION BETWEEN MOMENTS AND CUMULANTS 

Moments and cumulants which enters the expansion of exponent 

{oo k ^ oo k 

are related to each other as |21|] 

Zi = Fi 

Z2 = F2 -\- Fi 

= F^ + 3F1F2 + F^ 
Z^ = F^ + AF1F3 + 3F| + 6Ff F2 + F^ 



r=l 



ki\ J \kr\ J ill . . . i 



where summation is over all positive numbers and different positive numbers 

{fci, . . . , ky.} such that kiii + . . . + k^ir = k. 

APPENDIX C: ELLIPTIC THETA FUNCTIONS 

We adopt the following definition of the elliptic ^-functions: 

nez ^ ^ 

= ri{T) exp |7rir(a;^ — a + |) + 27rz(i — a){z + | — /?)} 

00 



00 

X 

n=0" 
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These should be compared with the notations of Mumford . 
The elhptic ^-functions satisfies the heat equation 



d_ 



a,(3 



1 92 



Relation to standard notations is 



(CI) 



dQfl{z,T) = 6*1 (2, r) 
Oo,^{z,t) = 92{z,r) 
Oi,o{z,t) = 9^{z,t) 
e^^^{z,T) = 9,iz,r) 



The Dedekind ?7-function is usually defined as 



00 



2niTn 



n=l 



Considering the functions Oa^pir) = 9a,i3{0,T) and ?7(r) of pure imaginary aspect ratio, 
T = ip, we obtain the identity 



In 



r]{ip) 



vrp = X] In I 1 - e-2-p("+")-2«/3 



n=0 



_^ ^ In I ]_ _ g-27rp(n+l-a)-27ri/3 



n=0 



(C2) 



Taking the derivative of Eq. ( P2|) with respect to p we can obtain the following useful 
identity 



00 00 



n=0 m=l 



[n + a)e 



-27rm[ p(n+«)+i/3 ] 



+ in + 1 — a)e 



-27rm[p(n+l-«)+i/3] 



— H 7^ In 

2 27: dp 



r]{ip) 



(C3) 



Taking the second derivative of In \9a,/3{z, ip)\ with respect to z at z = and using the heat 
equation Eq. (|Cl| ) we obtain 



^ ' - I -2nm[p{n+a)+il3\ 



'-^H \ n=Om=l 



1 



n + a 



-2-^^ln\e^A^p)\ + ^(-^ln\eaA^p)\ 



n + 1 — a 
\ 2 



-27rm[p(ra+l-a)+i/3] 



(C4) 
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APPENDIX D: KRONECKER'S DOUBLE SERIES 



Kronecker's double series can be defined as ||23[| 



p\ ^-2-Ki{na+m(3) 



{m,n)7^{0,0) 

In this form, liowever, tliey cannot be directly applied to our analysis. We need to cast 
them in a different form. To this end, let us separate from the double series a subseries with 
m = 

K"'/5(^) = Pi — y ^ Pi — y y ^ 

Here the first sum gives nothing but Fourier representation of Bernoulli polynomials Eg. ([A^) . 
The second sum can be rearranged with the help of the identity 

which can easily be derived from Eg . (|A4|) differentiating it p times. The final result of our 
resummation of double Kronecker sum is 

oo 

m^O n=Q 

Considering real part of the Kronecker sums with pure imaginary aspect ratio, r = ip, we 
can further rearrange this expression to get summation only over positive m > 1 



/o^Xl oo — 27ri(na+m/3) _i_ „— 27ri(?i(l— a)— m/3) 

Bl - ReK^^(zp) = Re V V , ^ 



oo oo 

a] 

Km=l n=0 

oo oo 



m=l n=0 



APPENDIX E: ASYMPTOTIC EXPANSION OF THE DIGAMMA FUNCTION 

Let us start with well known expansion of the digamma function ip{N) ||T 



^Ij{x) = Inx- 



^ B2p 1 



2x ~^ 2p x'^P 

oo R 1 

= lnx-E(-ir^^ (El) 

Plugging in the above expansion x = N + a and expand the resulting factors In (1 + a/N), 
(1 + a/N)~P in powers of N~'^ we obtain 
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ij(N + a) = iniv - y (-i)p^ - y y(-i)fe+PB /^^ , , 

PNV ^tl^o^ ' " k\p\ NP+' 

1=1 p=0 



\l-p)\p\ 

Using the relation between Bernoulli polynomials and Bernoulli numbers Bp 

l\ 



a 



i-p 



we finally obtain the Eq. (|2! 



OO go -1 

^(N + a) = lnN- Ti-lY^^ 



(E2) 



(E3) 



(E4) 



APPENDIX F: REDUCTION OF KRONECKER'S DOUBLE SERIES TO THETA 

FUNCTIONS 

Let us consider Laurent expansion of the Weierstrass function 



1 



+ E 



1 



TjZ' 



{z — n — rm)2 (n + rm)^ 

2p~2 



p=2 



The coefficients dpir) of the expansion can all be written in terms of the elliptic 6'-functions 



with the help of the recursion relation ^4 

3 



P (p-3)(2p+l) 
where first terms of the sequence are 



(a2ap_2 + ctsctp-s + • • • + 0^-202) 



02 - Y5 - 626I + 6pj) 
ag = ^(2a^ + 3a2) 



3*^2 
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Kronecker functions K^'p (r) are related directly to the coefficients ap{T) 

T^o,0/ X ^ (2p)! ap(r) 
(-47r2)P (2p - 1) 

Kronecker functions K2p'^(r) with a and /? taking values or 1/2 can in their turn be related 
to the function (r) by means of simple resummation of Kronecker 's double series 

Kj^(r) = 2K°'°(2r)-K;'°(r) 

K|'°(r) = 2i-^K;'0(r/2)-K°'0(r) 

K|'^(r) = (l + 22-f)KjO(r) -2^-fKjO(r/2) -2Kj'°(2r) 

Thus, Kronecker functions K2p^(r) with a and /3 taking values or 1/2 can all be expressed 
in terms of the elliptic ^-functions only. For practical calculations the following identities 
are also helpful 

2^i(2T) = el - el ei{T/2) = 2^2^3 

2ei{2T) = ei + el ei{T/2) = ei + ei 

2^|(2r) = 2^3^4 el{T/2) = el - el 

From the general formulas above we can easily write down all the Kronecker functions 
that have appeared in our asymptotic expansions 

Kr(r) = u^n - m - e\et) 

^'^{^) = i-M^e\){e\-e\){e\^et) 

K|'°(r) = -^(^l + ^l)(i^! + « 

n^{r) = ^,{ei-et){%el-e\et) 

Note that when p — cxd we have limits ^2 — > 0, ^4 — > 1, ^3 — > 1 and the Kronecker's function 
reduce to the Bernoulli polynomials. 
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